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Abstract: We provide a novel procedure to obtain complex PT-symmetric multi-particle 
Calogero systems. Instead of extending or deforming real Calogero systems, we explore 
here the possibilities for complex systems to arise from real nonlinear field equations. We 
exemplify this procedure for the Boussinesq equation and demonstrate how singularities in 
Q-ij real valued wave solutions can be interpreted as N complex particles scattering amongst 

r"| \ each other. We analyze this phenomenon in more detail for the two and three particle 

case. Particular attention is paid to the implcmcntion of PT-symmetry for the complex 
multi-particle systems. New complex PT-symmetric Calogero systems together with their 
• classical solutions are derived. 

o 

^ ! 1. Introduction 

ON . 

The analytic continuation of real physical systems into the complex plane is a principle 

which has turned out to be very fruitful, since many new features can be revealed in this 
manner which might otherwise be undetected. A famous and already classical example, 
proposed more than half a century ago, is for instance Heisenberg's programme of the 
analytic S-matrix PJ. Here our main concern will be complex multi-particle Calogero 
systems, in particular those exhibiting PT-symmetry 

Quantum systems are said to be PT-symmetric when they are invariant under si- 
multaneous parity P and time reversal T transformations. When the Hamiltonian, not 
necessarily Hermitian, exhibits this symmetry, i.e. [H,VT] = 0, and moreover when all 
wave-functions are also invariant under such an operation this property is referred to as 
unbroken PT-symmetry. The virtue of this feature is that it is a sufficient property to 
guarantee the spectrum of the Hamiltonian H to be real. The underlying mechanisms 
responsible for this are by now well understood f|[ f|, ||, |(| and may be formulated al- 
ternatively in terms of pseudo/quasi-Hermiticity; for definitions see for instance @ and 
references therein. 

There are two fundamentally different possibilities to view complex systems: One may 
either regard the complexified version just as a broader framework, as in the spirit of the 
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analytic S-matrix, and restrict to the real case in order to describe the underlying physics 
or alternatively one may try to give a direct physical meaning to the complex models. 

With the latter motivation in mind complex 'PT-symmetric Calogero systems have 
been introduced and studied recently || |9|, [ji], [ll], 12]. The hope for a direct physical 



interpretation stems from the fact that unbroken PT-symmetry will guarantee their eigen- 
spectra to be real and allows for a consistent quantum mechanical description, i.e. such 
systems constitute well defined quantum systems which have been overlooked up to now. 
Nonetheless, so far any such proposal lacks a direct physical meaning and the complexi- 
fications are generally introduced in a rather ad hoc manner. Here our main purpose is 
to demonstrate that various complex Calogero models appear rather naturally from real 
valued nonlinear field equations and thus we provide a well defined physical origin for these 
systems. 

The solutions for the real Calogero systems were found in the reverse order when 
compared to the usual way progress is made, i.e. the quantum theory was solved before the 
classical one. Calogero solved first the quantized one-dimensional three-body problem with 
pairwise inverse square interaction [|l3] and subsequently constructed the ground state of 
the iV-body generalization [14] described by the Hamiltonian 



N 2 - N 



i=l i^j 



with g G R being the coupling constant. Marchioro [15, [L6|] investigated thereafter the 



classical analogue of these models obtaining a solution to which we will appeal below. The 



integrability of these classical counterparts was established later by Moser [17], using a Lax 
pair consisting of matrices L, M, with entries 

Lij = Pi5ij + -^-(1 - %), (1-2) 

^k) y^i •bj) 

constructed in such a manner that the Lax equation 

§ + [M,L]=0 (1.4) 
becomes equivalent to the Calogero equations of motion, 



We use the notation i = y/—T throughout the manuscript and abbreviate time derivatives 
as usual by dxi/dt = and d 2 Xi/dt 2 = X{. Integrability follows in the standard fashion by 
noting that all quantities of the for I n = tr(L n ) jn are integrals of motion and conserved 
in time by construction. 
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Calogero systems have become very important in theoretical physics, having been 
explored in various contexts ranging from condensed matter physics to cosmology, e.g. 
fl8| , [lS| , p0| . The main focus of our interest here are the complex extensions which have 
been studied recently in connection with PT-symmetric models || ||, |H], [II], 12 1. 



The idea of exploiting PT-symmetry in order to obtain models with real energies can 
be adapted to classical systems as well and has been used to formulate various complex 
extensions of nonlinear wave equations, such as the Korteweg-de Vries (KdV) and Burgers 



equations |21, 22, 23, pj]. In the classical case the reality of the energy is ensured in an even 



simpler way, as in that case the PT-symmetry of the Hamiltonian is sufficient. Remarkably 
these systems allow for the existence of solitons and compacton solutions p5| , 26]. 

Here we shall explore PT-symmetry in a context where the complex extensions or 
deformations do not need to be imposed artificially, but instead we investigate whether 
this symmetry is already naturally present in the system, albeit hidden. To achieve this 
goal we exploit the fact that nonlinear equations, such as Benjamin-Ono and Boussinesq, 
can be associated to Calogero particle systems. We explore these connections and are then 
naturally led to complex PT-symmetric Calogero systems. 

In the next section we shall demonstrate how a complex one-dimensional PT-symmetric 
Calogero system is embedded in a real solitonic solution of the Benjamin-Ono wave equa- 
tion and how constrained PT-symmetric Calogero particles emerge from real solutions of 
the Boussinesq equation. Thereafter we construct the explicit solution of the three-particle 
configuration with the aforementioned constraint and show that the resulting motion, un- 
like in the unconstrained situation, cannot be restricted to the real line. We shall also 
establish that a subclass of this constrained Calogero motion is related to the poles in the 
solution of different nonlinear KdV-like differential equation. The relation of these com- 
plex particles with previously obtained PT-symmetric complex extensions of the Calogero 
model |^7| is discussed in section 5, where we demonstrate that they are different from 
those proposed here. Our conclusions are drawn in section 6. 



2. Poles of nonlinear waves as interacting particles 

The assumption of rational real valued functions as multi-soliton solutions of nonlinear 
wave equations was studied more than three decades ago by various authors, see e.g. |^8| . 
We take some of these findings as a setting for the problem at hand. In order to illustrate 
the key idea we present what is probably the simplest scenario in which corpuscular objects 
emerge as poles of nonlinear waves, namely in the Burgers equation 

u t + au xx + f3{u 2 ) x = 0. (2.1) 

Assuming that this equation admits rational solutions of the form 

. . la 1 , . 
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it is straightforward to see that surprisingly the N poles interact with each other through 
a Coulombic inverse square force 

N 1 

Xi(t) = — 2a y 7 — — tttt . (2.3) 

W ^[^(*)-^(*)] 2 ^ ' 

This pole structure survives even after making modifications in the ansatz for the wave 
equation, although the nature of the interaction may change. By acting on the second 
derivative in Burgers equation with a Hilbert transform 

i r +oc u(z) 

Huix) = -PV / dz^^, (2.4) 

7T y_ 00 Z-X 

we obtain the Benjamin-Ono equation 3C] 

u t + aHu xx + (3(u 2 ) x = 0. (2.5) 

As shown in [p]}l , the ansatz proposed for the equation above which will allow for 
similar conclusions has a slightly different form, 

N / \ 

being, however, still a real valued solution with the only restriction that the complex poles 
satisfy complex Calogero equations of motion 

N 1 



Note that there is a difference in the power laws appearing in ( [2.3; ) and (2/7), but 
more importantly that equation ( |2.2[) has real poles, whereas ( |2.6| ) has complex ones. We 
stress once more that the field u(x, t) is real in both cases. Hence, this viewpoint provides 
a nontrivial mechanism which leads to particle systems defined in the complex plane. 

Interesting observations of this kind can be made for other nonlinear equations as well, 
but not always will the ansatz work directly, that is without any further requirements as in 
the previous cases. In some situations additional conditions might be necessary. Examples 
of nonlinear integrable wave equations for which such type of constraints occur are the 
KdV and the Boussinesq equations, 

Ut + (au xx + Pu 2 ) x = and u tt + (au xx + (3u 2 - ju) xx = 0, (2.8) 

respectively. For both of these equations one can have "iV-soliton" solutions 1 of the form 

N 



a \ 1 

U[X,t) = -6— } 7 T-TTTT, (2.c 



1 Soliton is to be understood here in a very loose sense in analogy to the Painleve type ideology of 
indistructable poles. In the strict sense not all solution possess the iV-soliton solution characteristic, that 
is moving with a preserved shape and regaining it after scattering though each other. 
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as long as in each case two sets of constraints are satisfied 

N N 

x k {t) = -UaJ^Mt) ~ Xj{t)y 2 , = ^2(x k (t)-x J (t))- 3 , (2.10) 

jj^k j^k 

and 

N N 

Zfc(t) = -24a^(x fc (t) - Xjit))- 3 , ± k (t) 2 = \2aY,{xk{t)-x 3 {t))- 2 + 7, (2.11) 

j^k jjtk 

respectively. Naturally these constraints might be incompatible or admit no solution at all, 
in which case ( ^Sj) would of course not constitute a solution for the wave equations (^) . 
Notice that if the x k (t) are real or come in complex conjugate pairs the solution ( |2.9[) for 
the corresponding wave equations is still real. 

Airault, McKean and Moser provided a general criterium, which allows us to view these 
equations from an entirely different perspective, namely to regard them as constrained 
multi-particle systems p8| : 

Given a multi-particle Hamiltonian H(x\, ...,xn,&i, with flow X{ = dH/dxi 

and Xi = —dH/dxi together with conserved charges I n in involution with H , i.e. vanishing 
Poisson brackets {H, I n } = 0, then the locus of grad(I n ) = is invariant with respect to 
time evolution. Thus it is permitted to restrict the flow to that locus provided it is not 
empty. 

Taking the Hamiltonian to be the Calogero Hamiltonian He it is well known that one 
may construct the corresponding conserved quantities from the Calogero Lax operator (ft. .2D 
as mentioned from I n = tr(L n )/n. The first of these charges is just the total momentum, 
the next is the Hamltonian followed by non trivial ones 

N 1 N N 

l!=S> h = H c (g), /s^E^EtI^V--- (2-12) 

According to the above mentioned criterium we may therefore consider an 13-flow 
restricted to the locus defined by grad(^2) = or an /2-flow subject to the constraint 
grad(/3 — 7/1) = 0. Remarkably it turns out that the former viewpoint corresponds exactly 
to the set of equations ( 2. 10]) , whereas the latter to ( [2.11 ) when we identify the coupling 



constant as g = —12a. Thus the solutions of the Boussinesq equation are related to the 
constrained Calogero Hamiltonian flow, whereas the KdV soliton solutions arise from an 
/3-flow subject to constraining equations derived from the Calogero Hamiltonian. 

As our main focus is on the Calogero Hamiltonian flow and its possible complexi- 
fications we shall concentrate on possible solutions of the systems (2.11) and investigate 



whether these type of equations allow for nontrivial solutions or whether they are empty. It 
will be instructive to commence by looking first at the unconstrained system. The classical 
solutions of a two-particle Calogero problem are given by 



•n..2(/) - 2R(I) + \/| + 4£(t-t ) 2 , (2.13) 
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with E, to being initial conditions and R(t) = the centre of mass velocity. Relaxing this 
condition by allowing boosts will only shifts the energy scale since the total momentum 
is conserved. Depending therefore on the initial conditions we may have either real or 
complex solutions. 

The three particle model, i.e. taking N = 3 in ( |1.1| ), is slightly more complicated. 
Marchioro [15| found the general solution by expressing the dynamical variables in terms 
of Jacobi relative coordinates R, X, Y in polar form via the transformations R(t) = 
{x l (t)+x 2 {t)+x 3 (t))/3, X(t) =r{t)sm(f>{t) = { Xl {t)-x 2 {t))/\/2 and Y(t) = r(t) cos </>(t) = 
(x 1 (t)+x 2 (t)-2x 3 (t)) /V6. The variables may then be separated and the resulting equations 
are solved by 



xi, 2 (t) = R(t) + -pr(£) cos 4>{t) ± -pr(t) sin <f>(t), 
V6 V2 



x 3 (t) = R(t) 



r(t) cos <j)(t), 



(2.14) 
(2.15) 



where 



R(t) — Rq + tRQ 
r(t) = 



<Kt) 



B 2 



E 



- cos 
3 



+ 2E{t -t ) 2 , 



<f sm 



sin 1 (ip cos30 o ) — 3 tan 1 ( ^— — (t — to) 

V B 



(2.16) 
(2.17) 

(2.18) 



The solutions involve 7 free parameters: The total energy E, the angular momentum type 
constant of motion B, the integration constants to, (j> Q , Rq, Rq and the coupling constant 
g, with the abbreviation cp Q = \/l — 9g/2B 2 . We note that, depending on the choice of 
these parameters, both real and complex solutions are admissible, a feature which might 
not hold for the Calogero system restricted to an invariant submanifold. 

Let us now elaborate further on the connection between the field equations and the 
particle system and restrict the general solution ( 2.14j )-( p.l6| ) by switching on the additional 



constraints in (|2.11[) and subsequently study the effect on the soliton solutions of the 
nonlinear wave equation. Notice that the second constraint in ( [2.11 ) can be viewed as 
setting the difference between the kinetic and potential energy of each particle to a constant. 
Adding all of these equations we obtain He = N*y/2, which provides a direct interpretation 
of the constant 7 in the Boussinesq equation as being proportional to the total energy of 
the Calogero model. 



3. The motion of Boussinesq singularities 

The two particle system, i.e. N = 2, is evidently the simplest /2-Calogero flow constrained 
with grad(/3 — jli) = as specified in ( |2.11| ). The solution for this system was already 
provided in pSfl , 

x li2 (t) = K± y /'y(t-K) 2 -3a/'y, (3.1) 
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with k, k taken to be real constants. In fact this solution is not very different from the 
unconstrained motion shown in the previous section ( 2.1 3| ) . The restricted one may be 
obtained via an identification between the coupling constant and the parameter in the 
Boussinesq equation as k = 2R(t), E = 7/4, k = to and g = —3a/ 4. The two soliton 
solution for the Boussinesq equation (|2.9|) then acquires the form 



u(x, t) 



-12* 7 7(X 



k? 



3a 



p '[ 7 (x-K) 2 -7 2 (t-K) 2 + 3a] 2 ' 



(3.2) 



which, in principle, is still real-valued when keeping the constants to be real. When in- 
specting Q3.2] ) it is easy to see that the two singularities repel each other on the x-axis 
as time evolves, thus mimicking a repulsive scattering process. However, we may change 
the overall behaviour substantially when we allow the integration constants to be complex, 
such that the singularities become regularized. In that case we observe a typical soli- 
tonic scattering behaviour, i.e. two wave packets keeping their overall shape while evolving 
in time and when passing though each other regaining their shape when the scattering 
process is finished, albeit with complex amplitude. A special type of complexification oc- 
curs when we take the integration constants k, k to be purely imaginary, in which case 
Q3.2D becomes a solution for the PT-symmetrically constrained Boussinesq equation, with 
VT : x — ► —x,t — > —t,u — ► u. We depict the described behaviour in figure 1 for some 
special choices of the parameters. 
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Figure 1: Time evolution of the real part of the constraint Boussinesq two soliton solution (3.2) 
with k = «2.3 , k = — lO.Q, a = —1/6,(3 = 5/8 and 7= 1. 
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For larger numbers of particles the solutions have not been investigated and it is not 
even clear whether the locus of interest is empty or not. Let us therefore embark on solving 
this problem systematically. Unfortunately we can not simply imitate Marchioro's method 
of separating variables as the additional constraints will destroy this possibility. However, 
we notice that (2.11) can be represented in a different way more suited for our purposes. 



Differentiating the second set of equations in ( p. 11 ) and making use of the first one, we 
arrive at the set of expressions 

A (x k (t) + Xj(t)) 

which are therefore consistency equations of the other two. 

We now focus on the case N = 3. Inspired by the general solution of the unconstrained 
three particle solution ( 2. 14] ) and ( p. 15 ), we adopt an ansatz of the general form 



xi,2(t) = A (t) + A 1 (t) ± A 2 (t), (3.4) 
x 3 (t) = 4>(t) + AAi(t), (3.5) 

with Ai(t), i = 0,1,2 being some unknown functions and A a free constant parameter. 
We note that A ^ 1, since otherwise the three coordinates could be expressed in terms of 
only two linearly independent functions, Ao(t) + Ai(t) and A 2 (t), and we would not able 
to express the normal mode like functions Ai(t) in terms of the original coordinates Xi(t). 
Calogero's choice, A = —2, in equation ( 2.15[ ), allows an elegant map of Cartesian coordi- 



nates into Jacobi's relative coordinates, but other possibilities might be more convenient 
in the present situation. Here we keep A to be free for the time being. 

Substituting this ansatz for the Xi(t) into the second set of equations in ( |2.11| ) and 
using the compatibility equation (|3.3j), we are led to six coupled first order differential 
equations for the unknown functions Ao(t), Ax(t), A 2 (t) 

2g + 2A~W + 2AZW = ' ( } 

(A (t) + A 1 (t)±A 2 (t)) 2 - 1 , 1 1 



2g 8A 2 (t)i 2A T (t)? 

2Ao(t) + (A + l)ii(i) + A 2 (t) 2A Q (t) + (A + l)ix(t) - A 2 (t) 



+" VTTT^ + ^T^2 = 0, (3.7) 

0, (3.8) 



A-(t) 3 A+(tf 

Ap^ + A^t) 2A (t) + (\ + i)A 1 (t)±A 2 (t) 

^Mtr + a^w " ( } 

For convenience we made the identifications A±(t) = A 2 (t) ± (A — l)Ai(t). 

From the latter set of equations above, ([T^) and fl3.9|) , we can now eliminate two of 
the first derivatives together with the use of the conservation of momentum. Depending 
on the choice, the remaining Aj(t) are eliminated with the help of the first three equations 
J3.6| ) and ( |3.71 ). The two equations left then become multiples of each other depending only 
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on Ai(t) and ^(i). Subsequently we can express A2(t), and consequently Ao(t), Ai(t), in 
terms of Ai(t) as the only unknown quantity. In this manner we arrive at 



A 2 (t) 
A (t) 



V-5-4 7 (A-l)M 1 (t)2 



2^ 



3gy7(2 + A) 



( A _i)[ 5 + i6 7 ( A _ IfA^tfY 



(1-A)b + 16 7 (A-1)M 1 (^ 



(3.10) 
(3.11) 
(3.12) 



with g = —12a. This means that once we have solved the differential equation ( |3.12j ) for 
Ai(t) the complete solution is determined up to the integration of Ao(t) in ( [3.11| ) and a 
simple substitution in ( 3.1C| ). In other words we have reduced the problem to solve the set 
of coupled nonlinear equations ( |2 . 1 1| ) to solving one first order nonlinear equation. 

Let us now make a comment on the number of free parameters, that is integration 
constants, occurring in this solution. In the original formulation of the problem we have 
started with 3 second order differential equations, so that we expect to have 6 integration 
constants for the determination of x±,X2 and X3. However, together with the additional 3 
constraining equations this number is reduced to 3 free parameters. Finally we can invoke 
the conservation of total momentum from ( |2.12| ), which yields 3Ao(t) + (A + 2)A\{t) = 
and we are left with only 2 free parameters. We choose them here to be the two arbitrary 
constants attributed to the integration of ^4o(^) m (3.11) and Ai(t) in ( p. 12 ), respectively. 

In turn this also means that, without loss of generality, we may freely choose the 
constant A introduced in Q3,5| ). Indeed, keeping it generic we observe that the solutions for 
the Ai(t) do not depend on it despite its explicit presence in the equations ( 3.10| ), ( |3.11| ) 
and ( 3.12p . The most convenient choice is to take A = —2 as in that case the equations 
simplify considerably. 

Let us now solve ( 3.10| ), ( 3.11 ) and ( 3.12 ) and substitute the result into the original 
expressions (|3.4| ) and ( ^ ) in order to see how the particles behave. We find 



Xl ^ t) = C0 + ^ t + T2{w) 



xs(t) = C + y/jt - - 



g 



7 



4^ U(*) 



9 + m 



6 U(i) 7 

where for convenience we introduced the abbreviation 



^r{^/igt + ci) + J 9 3 7 3 + [54 7 2 (V7S* + ci)] 2 



(3.13) 
(3.14) 

(3.15) 



The above mentioned two freely choosable constants of integration are denoted by cq and c\. 
As in the two particle case, we may once again compare this solution with the unconstrained 
when considering the Jacobi relative coordinates 



one in fljTg), (gji 

R(t) = c + t^f, 



g_ 

67 



and tan(/>(i) 



: 91 + t\t) 
' 91 -fit)' 



(3.16) 



-9- 



From real fields to complex Calogero particles 



We observe that the solution is now constrained to a circle in the Xl"-plane with real 
radius when gj € Mr. The values for (f>(t) lead to the most dramatic consequence, namely 
that the particles are now forced to move in the complex plane, unlike as in unconstrained 
Calogero system or the N = 2 case where all options are open. 

Interestingly, despite the poles being complex, we may still have real wave solutions 
for the Boussinesq equation. Provided that 7, g, cq, c\ £ M the pole xs(t) is obviously 
real whereas x±(t) and X2(t) are complex conjugate to each other, such that the ansatz 
fl2.g| ) yields a real solution 

u(x,t) = — £ ^+ (3.17) 

"MGfe-f)) 

gW - 12<ry V(*) ~ 4 7 (18 7 y? 2 - 9)m 2 + 12 7 ^(t) 3 + £(i) 4 



216a 9j _. , 9 



( 5 V + 697V? (t) + + <j)«S) 2 - 6 7 y«*) 3 + ?(*) 4 ) 2 



with ip = Co + y / 7< — x. 

Due to the non-meromorphic form of £(t) it is not straightforward to determine how 
the solutions transforms under a 'PT-transformation. Nonetheless, the symmetry of the 
relevant combinations appearing in (ETIaD and ( |lg) can be analyzed well for cq,c\ & 



and 7 > . In that case the time reversal acts as T : i i ^|^y i > which 

implies "PT : — > —xi(t) for i = 1,2,3. Thus, the solutions to the constrained problem 
are not only complex, but in addition they can also be PT-symmetric for certain choices 
of the constants involved. 

4. Different types of constraints in nonlinear wave equations 

It is clear from the above that the class of complex (PT-symmetric) multi-particle systems 
which might arise from nonlinear wave equations could be much larger. We shall demon- 



strate this by investigating one further simple example which was previously studied in [33] 
and also refer to the literature [|32| for additional examples. One very easy nonlinear wave 
equations which, because of its simplicity, serves as a very instructive toy model is 

u t + u x + u 2 = 0. (4.1) 

We may now proceed as above and seek for a suitable ansatz to solve this equation, possibly 
leading to some constraining equations in form a multi-particle systems. Making therefore 
a similar ansatz for u(x,t) as in ( |2.6| ) or (|2.9| ) we take 



It is then easy to verify that this solves the nonlinear equation ( ^4.1| ) provided the Zi(t) obey 
the constraints 

~f- Zi{t)-Zj{t) 



- 10 - 



From real fields to complex Calogero particles 



We could now proceed as in the previous section and try to solve this differential equation, 
but in this case we may appeal to the general solution already provided in 33], where it 
was found that 

l + tf(x-t) 

solves (|4.1|) for any arbitrary function f(x) with initial condition u(x, 0) = f(x). Com- 



paring ([0|) and (4.2) it is clear that the Zi(t) can be interpreted as the poles in (p~4|), 
which becomes singular when x — > Zj (t) = t + /^(-l/t), with % G {1, N} labeling the dif- 
ferent branches which could result when assuming that / is invertible but not necessarily 
injectively. Making now the concrete choice for / to be rational of the form 



N 



f( x ) = ^2 — > witl1 a i, a i G c > ( 4 -5) 

1=1 a i~ X 

we can determine the poles concretely by inverting this function. First of all we obtain 
from the initial condition that 

Zi(0)=ai, Zi(0) = l + ai and %(0) = V 0< ° i . (4.6) 

*r?.a.i - ay 



The first two conditions simply follow from the comparison of ( f4.4| ) and ( |4.2| ), but also 
follow, as so does the latter, from taking the appropriate limit in ([4.51). We note that the 



total momentum is conserved for this system Y2i=i ^»(*) = N + Y^i=i a i- 

Let us now see how to obtain explicit expressions for the poles. Inverting (4.5) for 
N = 2 it is easy to find that for generic values of t the poles take on the form 



zi, 2 (t) = t + ^- + ^j-t±-^/al 2 + 2a l2 a l2 t + a 2 12 t 2 , (4.7) 

where we introduced the notation ajj = a?j — ay, otij = ai + aj and analogously for a — > a. 
We note that in the case N = 2 the constraint (|4.3| ) can be changed into two-particle 
Calogero systems constraint with the identification g = a\a 2 a\ 2 - 

Next we consider the case N = 3 for which we obtain the solution 

zi{t) = t-^f- + s+(t) + 8-(t) (4.8) 
Z2fl® = t - ^ - \ Wit) + ._(*)] ± [s + (f) - «_(*)] , (4.9) 



where we abbreviated 



1/3 



s±(t) = [r(t) ± vV 2 (*)+<Z 3 (*)J , (4-10) 

r(t) = 9a(t)^)-27c(t)-2a 3 (t) ^ ^ = 3fr(t) - a 2 (t) ^ ^ 

a(t) = -at — a 2 — a 3 — t{a\ + a 2 + a 3 ), (4.12) 

b(t) = a\a 2 + «2«3 + «i«3 + i[aia23 + a2«3i + 03021], (4-13) 

c(t) = —t{aia 2 a2, + a 2 cn2,Oii + aj,a\a 2 ) — aia 2 a^,. (4-14) 
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In terms of Jacobi's relative coordinates this becomes 

R{t) = t-\a{t), r 2 (t) = 6s + (t)s.(t) and tanflt) = i *"® ~ S+ ® , (4.15) 
3 s-(t) + s+(t) 

which makes a direct comparison with the constrained Calogero system ( |3.16j ) straight- 



forward. As the system ([4.8|), (4.9) involves more free parameters than the constrained 



Calogero system ( p. 16 ), we expect to observe some relations between the parameters Oj, en 



to produce the right number of free parameters. Indeed, we find that for 

^ = -^H(ai - aj )- 2 (4.16) 



2 

3T* 



and the additional constraints 

1 3 2 3 

co = -^2ai, ci = — Y[ (aj +a k -2ai), g = A^a i a j -a 2 , 7 = 1. (4.17) 

8=1 l<J<fc<3 1=1 

j,k^l i<3 

the two systems become identical. Thus we have obtained an identical singularity structure 
for two quite different nonlinear wave equations. 

5. Complex Calogero models and "PT-deformations 

We have demonstrated in section || that the solutions of the constrained Calogero models 
are intrinsically of a complex nature. As there have been various proposals before in the 
literature suggesting complex Calogero systems in form of 'PT-symmetric deformations, 
we shall now compare them with the above outcome. We will argue that the deformations 
presented here are new and different to those suggested up to now. 

The simplest "PT-deformation of any model is obtained just by adding a VT- invariant 
term to the original Hamiltonian. For a many-body situation, this was proposed for the 
first time in the framework of A n Calogero models by introducing the Hamiltonian || 

N 



H(q,p) = H c (q,p) + £ f ■ (5.1) 

(Xi Xj) 



In [10 1 it was shown that this simply corresponds to shifting the momenta in the standard 
Calogero Hamiltonian together with a re-definition of the coupling constant, which means 
the above construction is certainly quite different from the proposal (5.1). 



The second type of deformation [jE^] consists of replacing directly the set of ^-dynamical 
variables q = {q%, . . . ,qg} and their conjugate momenta p = {p±, . . . ,pi} by means of a 
deformation map e : (q,p) — > (q,p), whereby the map is constructed in such a way that 
the original invariance under the Weyl group W is replaced by an invariance under a VT- 
symmetrically deformed version of the Weyl group yV^ T . In terms of roots the map is 
defined by replacing each root a by a deformed counterpart a as e : a — > a, whereby 
the precise form of the deformation ensures the invariance under as specified in 
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fl2f |. Expanding the momenta as p = Y^i^iai, with G R, this means for the Calogero 
Hamiltonian 

e : H c {q,p) -> H VT (q,p) = ^^KiKjaidtj + ^ ^ g . (5.2) 

In order to find the concrete forms for q and p we need to be more specific about the 
algebras involved. Let us therefore examine the models based on the rank 2 algebras A2, 
B2 and G2. Depending on the dimensionality of the representation for the simple roots, 
we obtain either a two or a three particle systems and may therefore compare with the 
solutions found in the previous sections. In all cases the deformations of the simple roots 
ai and «2 take on the general form 

0.1(e) = R(e)a\ + H(e)Ki2\2, and a 2 (e) = R(e)ct2 — H(£)K 2 i\i, (5.3) 

with Ai, A2 being fundamental weights obeying 2A« • aj/a 2 = 5ij, the functions R(e), 1(e) 
satisfy lim e ^o R{e) = 1, hm e _>o = and Ky = 2aj • aj/a 2 are the entries of the Cartan 
matrix. Let us now take the following two dimensional representations for the simple roots 
and fundamental weights 



(5.4) 



We easily verify that this reproduces the correct entries for the Cartan matrices A 2 : 
K X i = K 22 = 2, K l2 = K 2 i = -1, B 2 : K n = K 22 = 2, K 12 /2 = K 21 = -1 and 
G2 '■ Ku = K 2 2 = 2, K\ 2 = K 2 i/3 = —1. Having constructed the deformed roots we 
compute next the deformed conjugate momenta and coordinates. In the representations 
fl5.4j ) the kinetic energy term changes just by an overall factor as 

p 2 = [R(e) - Vgl(e)] P 2 with va 2 = ub 2 = 1> V G 2 = — V3. (5.5) 



A 2 


a\ 


= (1, 


-V3), 


OL2 = 


(1,V3), 


Ai 


= |ai + |a2, 


A 2 


= \ct\ + |«2 


B 2 




= (1, 


-1), 


C*2 = 


(0,1), 


Ai 


= OL\ + <»2, 


A 2 


= \(X\ + «2, 


G 2 


a 1 


= (" 


V 2' V 1)1 


Ct2 = 


(V6,0), 


Ai 


= 3ai + 02, 


A 2 


= 3ai + 2a2- 



The specific choice R(e) = coshe and 1(e) = z^~ 2 sinhe, used in [12], keeps the kinetic 
energy term completely invariant, in the sense that the original and deformed momenta 
are identical. The dual canonical coordinates q are computed from 

a-q = q-a, a, q G R, a, q £ R iR. (5.6) 

We find 

qi = R(e)qi + iv g I(e)q2, and q 2 = R(e)q2 - w s I(e)qi. (5.7) 

We will now argue that ( |5.7| ) is always different from the constrained two particle solution 
of the Calogero model ( |3.1| ). In order to see this we recall first of all that for the solution 
to be 7'T-symmetric we require k, k £ iR. Equating now the sums x\ + X2 = q\ + q~2 we 
conclude that qi(t) = —q 2 (t) = —k/v s I(e) = const. Next we compute (x\ — X2) 2 , which 
yields 

j(t - R) 2 - 3q/7 = 2R 2 (e)q\(t). (5.8) 
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This equation is inconsistent as the right had side is real and time independent, a condition 
which can not be achieved for the left hand side. This proves our statement that the 



deformation method suggested here is genuinely different from the proposal in [12] in the 
two particle case. 

Keeping the deformed roots to be of the form ( |5.3| ), the three dimensional repesenta- 
tions for the simple roots 

A 2 :ai = (1,-1,0), a 2 = (0,1,-1), G 2 : e*i = (1,-1,0), a 2 = (-2,1,1), (5.9) 



yield the same result for the kinetic energy term ( |5.5|) , but obviously have to produce 
different dual canonical coordinates q. In this case we obtain 

q x = R(e) qi + i( s I(e)(q 2 - <fe), (5.10) 

q 2 = R(e)q 2 + i( g I(e)(q 3 - gi), (5.11) 

q 3 = R(e)q 3 + iC g /(e)(«i - 92), (5.12) 



where ( A2 = 1/3 and C Ga = — 1. Equating these solutions with ( |3.13| ), ( |3.14 ) and solving 



the resulting equations for the % with i = 1, 2, 3, it is easy to argue that the qi can not be 
made real, which establishes the claim that the solutions are also intrisically different for 
the three particle model. 



6. Conclusions 

Hitherto there have been two different types of procedures to complexify Calogero models. 
As explained in section 5 one may either add "PT-symmetric terms to the original Hamilto- 



nian which have turned out to be simple shifts in the momenta [1C] or one may directly 
deform the root system on which the formulation of the model is based [12]. In all these 
approaches the deformation is introduced in a rather ad hoc fashion. In this paper we have 
provided a novel mechanism, which has real solutions of physically motivated nonlinear 
wave equations as the starting point. The constrained motion of some solitonic solutions 
of these models then led to complex Calogero models, some of them being 'PT-symmetric. 

There are some obvious open problems left. For instance it would naturally be very 
interesting to study systems involving larger numbers of particles, which would correspond 
to higher soliton solutions for the nonlinear wave equations. Clearly the study of differ- 
ent types of wave equations, such as the KdV etc and their PT-symmetrically deformed 
versions would complete the understanding. 

Our analysis is schematically summarized in figure 2. 
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real multi-particle systems 




complex multi-particle systems 



complex nonlinear field equations 
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